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Ph' Abstract 

f^ i We present a systematic analysis of next-to-next-to-leading-order diagrams 

I that contribute to the pp — > ppvr" production at threshold. Analytic ex- 

C^ ' pressions are given for the effective transition operators, and the relative 

QQ ■ importance of various types of diagrams is discussed. The vertex-correction- 

TJ" . type graphs are found to give only small corrections to lower order graphs in 

^:;^ ! conformity with expectations. By contrast, we find very large contributions 

f^ I from two-pion exchange graphs that can be interpreted as a part of effective 

O^ ' (T- meson exchange diagrams. The recoil correction to the pion rescattering 
diagram also turns out to be large. 
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a. Introduction High-precision measurements |T|J^ of neutral pion production in proton- 
proton collisions pp — >■ ppir^ just above the threshold have generated renewed theoretical 
scrutiny of this reaction 0- [|T6|. This reaction is unique among two-nucleon pion produc- 
tion processes in that it is not well described by the single-nucleon process (Born term), 
Fig. 1(a), and the s-wave pion rescattering process. Fig. 1(b). The reason is that the "large" 
Weinberg- Tomozawa term does not contribute to the pp —>■ ppn^ reaction, in contrast with 
e.g. charged-pion production pp —^ pnTr~^, thus rendering pp —>■ ppir^ particularly sensitive to 
and hence an interesting testing ground for the less- well-understood "small" isoscalar s-wave 
pion rescattering terms. From the calculations done thus far it is clear that the two most 
basic processes. Fig. 1(a) and 1(b), give much smaller pp —>■ ppn^ cross sections than the 
measured values. Lee and Riska's model calculation [|] suggests that shorter-range isoscalar 
meson-exchange processes, like a- and u- exchanges between the two protons, might be very 
important in this reaction. 

In heavy-baryon chiral perturbation theory [HB^PT] [jT^ one can define the "large" and 
"small" terms in the Hamiltonian by way of chiral-order counting. In this language the first 
pion rescattering contributions to pp -^ ppn^ are one chiral-order higher than those of the 
charged-pion production (e.g., pp —>■ prni^). This means that pp —>■ ppn^ production is very 
sensitive to "correction terms", which tend to be "masked" by the leading terms in most 
low-energy pion-nucleon processes. Since the Born term and the pion rescattering term do 
not explain the pp — > pp-K^ cross section, we extend here our previous calculation fl^ to the 
next chiral order. 

It must be mentioned, however, that application of HB^PT to the A^A^ -^ NNtt pro- 
cesses is a delicate matter in at least two aspects. First, the non-negligible energy-momentum 
transfers involved can make less clear the distinction between reducible and irreducible dia- 



grams in Weinberg's chiral counting scheme. This has led Cohen et al. WM to propose new 
counting rules to be used for inelastic processes like A^A^ — > NNtt. It is noteworthy that 
in this modified scheme loop diagrams can be of the lowest order. Secondly, as the nucleon 
recoil involved becomes appreciable, the static nucleon "propagator", which is one of the 
basic features of the HB^PT formalism, can become increasingly problematic. Meanwhile, 
"improving" the HB^PT propagator by including a nucleon recoil term requires an extension 
of HBxPT from its original form. Despite the importance of these issues, we postpone the 
discussion of these matters to a future publication. In this Letter we rather use the standard 



counting rules a la Weinberg [|T8| within the framework of HB^PT. We concentrate on calcu- 
lations of the effective transition operators strictly within the narrow definition of HB^PT. 
We also relegate to the future the consideration of the initial- and final-state interactions 
in the transition amplitude, even though we are aware of the paramount importance of a 
full distorted- wave (DW) calculation for comparison with the experiment. This is, for one 
thing, because the transition operators resulting from our systematic treatment of the new 
diagrams exhibit very complicated energy- and momentum-dependencies, rendering a full 
DW calculation a highly non-trivial numerical task. Secondly, we wish to separate out the 
effects of higher-order diagrams in the "kernel" (irreducible part) of the reaction from the 
perhaps more mundane yet numerically important initial- and final-state interaction effects. 
The principal purpose of this note is to report analytic results for the transition operators 
that arise from a systematic treatment of the next chiral order diagrams in HB^PT. In order 
to gain some insight into the relative numerical significance of these diagrams, we compare 



their absolute values, at the threshold, to that of the one-pion exchange rescattering graph, 
Fig. 1(b). Our most important finding is that some of the two-pion exchange diagrams, 
which have no lower-order counterparts, give by far the dominant transition operators at 
the threshold (sometimes by an order of magnitude larger than the nominally "leading" 
ones). These new diagrams may be interpreted as a part of an effective a-meson exchange 
(see below). 

b. Conventions and other preliminaries The effective Lagrangian Cch in HB^PT is 
expanded as 



£,h = C'-^'> + C^''> + C^'^ + 



(1) 



£M represents a term of chiral order u with u = d + (n/2) — 2, where n is the number of 
fermion fields in the term, and d is the number of derivatives or powers of m.^. The explicit 
forms for the z/ = and 1 terms are lEDI : 
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where we have retained only terms of direct relevance to our present calculation. The 
SU(2) field matrix U{x) is non-linearly related to the pion field and has standard chiral 
transformation properties. We use the representation U{x) = Jl — [7t{x) / fj,]'^ + iT-7r{x) / f.^- 
as in Ref. [Q. N{x) represents the large component of the heavy-nucleon field; the four- 

:i, 0,0,0); D^N = [d^ + l[e,d^^])N isthe 



velocity parameter f^ is chosen to be v^ 



covariant derivative of A^; S^ is the (covariant) spin operator, which in the nucleon rest 

frame becomes S^ = (0,(t/2), and u^ = i[^^5^^ — ^d^j,^'^], where ^ = JU{x) [2^]. The low- 
energy constants Ci,C2 and C3 have been determined from other processes, see e.g. Refs. 
| r5| , |20[] . An explicit expression for C^'^\ which includes 0{m]^) recoil terms as well as terms 
containing new low-energy constants, can be found, e.g., in Ref. 



21 



In Weinberg's chiral counting |T^ each irreducible Feynman diagram carries a chiral 
order index v defined hy u = A — E^ — 2C + 2L + Y.i^ii where E^ is the number of nucleons 
in the Feynman diagram, L the number of loops, C the number of disconnected parts of the 
diagram, and the sum runs over all the vertices in the Feynman graph |T8|. In this note we 
consider irreducible diagrams with chiral orders up to z/ = 2 that give rise to pp —>■ ppTc^ 
transition operators. These diagrams are shown in Figs. 1-5. In the following T^'^' stands 
for a transition operator of chiral order u. The lowest-order transition operator for the 
Born diagram [Fig. (la)] has z/ = —1, and that for the rescattering diagram has z/ = +1 



[Fig. (lb)]. We define the first HB^PT calculation |1T3| , P^ which includes the aforementioned 
terms as the next-to-leading (NLO) calculation. Hence we decree our next chiral order [u = 
2) calculation to be next-to- next-to-leading (NNLO) order. With the use of Cch in Eq.(0), 
the two transition operators that feature in the NLO calculation are given (in momentum 
space) by ||T3| , [l^ 



where Pj and p- {j = 1,2) denote the initial and final momenta of the j-th proton. The 
four-momentum of the exchanged pion is defined by the nucleon four-momenta at the ttNN 
vertex: kj = pj — Pj, where pj = {Ep^,pj),Pj = {Epi,p-) with the definition Ep = (p^ + 

m'j^Y^'^ — rriN. The s-wave rescattering vertex function K{k,q) is calculated from Eq.(|^): 



^n ( r. I jA^^qh q-k 

-— 2Ci - C2 - ^ - C3^ 



.(*.,).^ 2c.-fe--^)^-c,i-f . (6) 



where k = {ko, k) and q = {uq, q) represent the four-momenta of the exchanged and final 
pions, respectively, and ujq = (g ^ + mlY^'^. 

c. Analytic results When we consider NNLO, we find 19 topologically distinct new types 
of diagrams that can potentially contribute to NN — > NNtt reactions. For a particular case 
of the pp -^ ppn^ reaction near threshold, the isospin selection rules and the s-wave character 
of the outgoing pion reduce this number from 19 to 6. We refer to these six as Type I, Type 
II, . . ., Type VI, and illustrate them in Figs. 2-4. Types I, II and III appear in Fig. 2, Type 
IV in Fig. 3, and Types V and VI in Fig. 4. Some of these diagrams have been considered by 
Gedalin et al. |^2[|, as we shall discuss later. We denote by Mi, Mn, ..., the NNLO transition 
operators that arise from the diagrams of Type I, II,..., respectively. In addition to these six 
operators, as will be explained in more detail below, there is an NNLO contribution from 
Fig. 5, which looks the same as Fig. 1(b), but which constitutes a higher-order correction 
to T*^"*"^) [eq.(||)]. We denote by Myn the transition operator due to this correction. Then 
the total NNLO transition operator we consider is given by 

T(+2) = Mi + Mil + Mm + Miv + My + Myi + Myn. (7) 

The explicit expressions for these operators are as follows. (In the expressions below the 
subscripts i and j refer to nucleon number 1 and 2.) 

^1 = 11 E ^.-fc.lXi + S (8) 






where 
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Xi = v-{pi+p'i) h{kj) +f-(g+^) v-{q+pi+p[) Ioiv-pi,-v-kj,k'^) + 

+ v-{q-^) v-{q-pi-Pi) Io{v-p[,vkj,k'^) (9) 

p. 
X2 = -v-{pi+p[) v-{5q-ki) v-kj I^{kj) - v-{q+^)v-{q+pi+p[) x 

X [Jo{v-pi)-jQ{v-{pi + kj)) + vkjV-{kj + 2pi)lQ{vpu-v-kj,k'^)] + 

-v-{q-^)v-{q-pi-Pi)[Jo{v-p'i)-Joiv-{p[-kj))+v-kjV-{kj-2p-)Io{v-Pi,v-kj,k'^j)]. (10) 
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where 



+ {-2ml + 2vp, v-p] + k]) I^{v.p'.,^.kj, e.)}. (11) 

^ni = |^ E [Si-kjX^ + 2te^,af3v''S^Stk''^X2], (12) 



Xl = -V-{pi+Pj+p'i+p'j)-In{kj) + 

H — / ._ ,A v-{Pi + (l+kj)yiv-p„-kj)-v-{p'j + q+kj)y{v-p'j,^j)] + 

^v-{Pi-q-kj)y{vp'i,kj)-v{pj-q-kj)y{vpj,kj)] (13) 



V 



{Pi-Pj] 



X2 = _ [v-{pi + q+kj)A{v-pi,-kj)-v{p'j + q+kj)A{v-p'j,-kj)] + 

'v-{p'i-q-kj)A{v-p[,kj)-v-{pj-q-kj)A{v-pj,kj)] . (14) 



v-iPr-Pj 
In Eqs.(|l3D and (|l|): 

y{uj,p) = 

^[4{vP)UP)~Muj)-Muj-vP) + {2ml-2{uj-vPf + P^)Io{uj,vP,P')] + 

^ n2ml-uj'-{uj-vPm ^^ 

\ 4P2 J 

x[2{vP)h{P)-Muj) + Jo{uJ-v-P) + {vP){2uj~vP)Io{uj,vP,P% (15) 

and 

4p2 

-ij„(^-^t,.P)+('^^^^^|^^) (J„M-J„(u;-«-P)) + 

x{(-^' + (^-P*)') Mv-Pi) + {'-ml + {vp'f) Mv-p',) + 
+v{p,+p[)A^ + [3kf-ml-2krq]X}, (17) 



where 
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h{ki) + 



+ m: 



il-v-piV-p^--kl\ Io{v-pi,v-ki,ki)--[Jo{v-pi) + Jo{v-p'i)]. 



(18) 
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x{v{pi+Pi)A^ + {{vpiy-ml)Jo{vp^) + {{v-p'i}'^-ml)Jo{vp'i)} 



(19) 



M, 



9A 



VI 



Sj ■ kj 



x{^^■(p^+pDA.+^;■(2g+pO t;-(-A;i + 2g+pO Jo{v-{p'i+q)) + 
+f ■(2g-pi) v-(-ki + 2q-pi) Mv ■ (pi- q))} . 



(20) 



As mentioned earlier, the graph in Fig. 5 with the relevant lowest order vertices gives T^"*"^) 
[eq.(|^)]. If we, however, use for the pion rescattering vertex in Fig.5 the lagrangian C^'^'^ ||21|| , 
there results an NNLO transition operator, Myn, which represents "recoil" term corrections 
of 0(m^^) to T(+^). Its explicit expression is 



M. 
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yA \ v^ //, N ^j'kjTj 
-r] 1^ i^i{kj,q) j2 T^— 



(21) 



where the expression for K,i{kj, q) is obtainable from Eq.(C.3) of Ref. |21|. 

The above expressions for the z/ = 2 diagrams contain four independent one-loop integrals, 
A^, Jo{uj), In{P^) and Io{uj, v-P, P^), of which the first three contain divergences. The finite 
parts are defined to include ln(m^/A). We choose the cut-off parameter to be A = 1 GeV. 
The two integrals, A^^ and Jo{uj), can be found in Ref. 
Feynman integral: 



while I-n^P"^) is a standard 
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The last integral, which is new, is given by 
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where s = m^-u;^ + x{2ujPo - P^) + x'^{P^ - P^), Pq = v-P and ^ = cu - xPq. Io{uj, v-P, P^) 



reduces to the integral 7o(f -P) given in appendix B of Ref. pO[ in the limit tu = and P^ 
= 0. This limit, however, is not applicable to the pp —>■ ppn^ reaction. 



Our NNLO diagrams contain the usual divergences which need to be regularized and 
then renormalized by appropriate counter terms. The single-nucleon process (Born term) 
of Fig.l(a) receives two types of corrections: (i) the finite "recoil corrections" of 0{m]^^) 
involving finite known parameters, and (ii) the (infinite) loop corrections. The u = 1 loop 
corrections, Tj'^.^\ to the Born amplitude were discussed in Ref. [|I^. We denote by T^^^^^^ the 
z/ = 2 loop and recoil corrections to the Born term. The recoil corrections to the Born term 
are reduced by (m^/mAr)'^''' . The divergences contained in the graphs in Fig. 2 are canceled 
( "renormalized" ) by counter-terms in C^'^^ corresponding to the five-point (uNNNN) vertex 
diagram (Fig.6). The same terms renormalize a part of the singularities in Mjv, Eq.(p!7|), 
coming from Fig.3. The remaining singularities in Miy are similar to those in the graphs in 
Fig. 4. To eliminate the latter singularities, £*^^^ must contain further counter-terms of the 



pion-nucleon scattering vertex type |^. This can be accomplished with the use of these 
counter-terms in graphs similar to the one in Fig. 5. We let T^^'^^^ stand for the v = 2 
transition operators that originate from such v = 2 counter-terms. The complete set of the 
z/ = 2 transition operators includes T^^l\ T^^r? ■, and T^o^^\ but we defer detailed discussion 



of these terms to a forthcoming paper |T^ and concentrate here on the finite parts of the 
following effective operators 

r = r^csc_^^(+2)_ (^24) 

d. Numerical results and discussion The purpose of this section is purely illustrative: we 
wish to have some idea as to the size of these corrections. A proper treatment of the derived 
transition operators must involve DW analyses, which we postpone to the forthcoming paper 
]TP| , where DW modifications to our numerical estimates described below as well as changes 
due to the use of a non-static nucleon propagator will be discussed. 

To proceed, we must fix the free parameters of the previous expressions. To NLO, as 
discussed in Ref. [0,|1^], the three parameters, Ci, C2 and C3 of Eq.(^, enter into the pion 



rescattering operator 7"^^**^. We shall use the three sets of parameters employed in Ref. [jTS 



Sets A, B and C in Table 1 summarize these values. Set A represents the central values of 
ci, C2 and C3 determined in Ref. pOf using the experimental values of the pion-nucleon a 



term, the nucleon axial polarizability a^ and the isospin-even s-wave ttN scattering length 
a"*". Sets B and C represent typical ranges of allowed values in the current determinations 
of theses parameters (see Ref. [Q for details). 

For simplicity we limit our consideration to the threshold kinematics, which means q^ = 
(r?7,7r,0) and the single exchanged boson (pion) of Fig. 1(b) has the four-momentum k'^ = 
(r7i^/2, k) with k"^ = — m^m^v- Then the final nucleon three momenta are p- = (for nucleon 
i = 1, 2), and n{k,q) in Eq.(Rp is fixed at |T3|P 



..,. ^(2,. _!(,, + , ,_j;d_)j. (25) 

Quantities of interest here are the magnitudes of the finite parts of the u = 2 transition 
operators, Eqs.(^-(^Ol), relative to the magnitude of the pion rescattering operator, 7"^'^*'^, 



^ Since we do not consider DW we leave out T^°^"' 



Eq. (D. Let us denote these ratios ^hj Rk = Mk/T^"''' (K = I, II, III, . . ., VII). Table 
2 gives Rk's for each of the parameter sets A, B, and C. We note that Rj, Rji and Rni, 
corresponding to the graphs in Fig. 2, give quite substantial individual contributions but 
Rji and Rjjj cancel each other at threshold.^ Most remarkably, Rjv corresponding to the 
pion-pion rescattering diagram. Fig. 3, is large, ranging 5 ~ 10. 

The appearance of these large individual contributions calls for an explanation. The 
two-pion exchange diagrams in Figs. 2 and 3 all involve one-loop integrals with three or four 
propagators. These loop integrals, for which typical four-momentum transfers k are large, 
can produce the factor k"^ = —m,T^m]\f in the numerator multiplying the integral I(){uj,v-P, P^) 
of eq. (^) . This factor turns out to be accompanied by some negative powers of fl coming 
from the vertices, resulting in a large enhancement factor, 0{k'^ f'"^ = —mT^mNf^'^ ~ —15). 
This feature essentially explains the large size of the two-pion exchange diagrams, although 
there are other diagram-specific numerical factors. 

The two-pion exchange diagrams in Fig. 2 can perhaps be viewed as a part of an effective 
cr-meson exchange that Lee and Riska [ffl found to be important in pp -^ ppix^ . It has been 



shown via soft-pion arguments [23| that the effective a- meson exchange can be understood 



as a two-pion exchange. The results in Ref. [23| lead us to suspect that the few diagrams 



we consider here are insufficient to generate the full strength of isoscalar two-pion exchange 
between two nucleons, but our NNLO results are indicative of the importance of the two- 
pion exchange diagrams for pp —>■ ppn^. Similar HBxPT two-pion exchange diagrams have 
been considered in calculating the scattering amplitudes for higher partial waves in A^A^ 



collision p4| . It is to be noted, however, that the higher partial wave amplitudes, which are 
only sensitive to peripheral NN scattering, can probe two-pion exchange contributions only 
for low three-momentum transfers. By contrast, in the A^A^ — >■ NNir reaction, the two-pion 
exchange diagrams are probed in a very different kinematical regime of "high" energy- and 
three-momentum transfers between the two nucleons. It is therefore not surprising that the 
roles of the two-pion exchange diagrams in our calculation are very different from those 



discussed in Ref. ||2^ 



The diagrams shown in Fig. 4 generate effective form factors at the pion-nucleon rescat- 
tering vertex in T^^^\ According to Table 2, the contributions of the corresponding opera- 
tors, Rv and Ryi, are less than 20% and 15% , respectively. Thus, the higher chiral-order 
vertex corrections to T^"*"^^ are found to be small, as expected from the general tenets of 
xPT. Meanwhile, Table 2 shows that Ryii = 1-9 ~ 2.6. This means that the combined 
pion recattering term is given by the sum of Fig. 1(b) and the Type VII contribution. 



Gedalin et al. p2| considered some NNLO diagrams within HB^PT. Numerically, they 
have found that the sum of Miv and My (in our notation) is large. This feature is confirmed 
by our numerical result. It is worth emphasizing, however, that of these two operators Miv 
is by far the predominant one. Although Gedalin et al. also considered Mn (our notation). 



^ This type of comparison is possible because, at threshold, effectively only one kind of spin 
operator appears in the pp -^ ppn^ transition operator. 

^ With DW this cancellation will not occur due to the different energy and three-monientum 
dependences. 



they left out Mi and Mm. According to Table 2, Mji and Mm are of equal importance, and 
their individual contributions are comparable to that of Miy. Our final remark is that our 
results contain a new loop-integral Iq{uj,v-p,p'^), Eq. p^), which does not appear in Ref. 

H- 

e. Conclusions We have evaluated the pp -^ ppir^ transition operators to NNLO in 
chiral expansion. It was found that xPT vertex corrections to the lower chiral order transi- 
tion operators, 7"^°™ and T^'^^'^ are indeed small. Thus, for this limited type of diagrams, 
xPT expansion seems to be under control. Meanwhile, the two-pion exchange contributions 
are found to be very large in our estimates. This result is consistent with the expectation 
that the pp -^ pprc^ reaction is sensitive to "heavy" -meson exchanges between nucleons. The 
phenomenologically important a-meson contributions @] seem to have discernible "represen- 
tatives" among the NNLO chiral perturbation diagrams considered here. It is not obvious 
whether one can interpret the large contributions from the individual graphs in Figs. 2 and 
3 as evidence for the non-convergence of the xPT expansion. These types of graphs make 
their first appearance only in the NNLO calculations, and therefore the convergence question 
can only be settled by calculating corrections to these NNLO diagrams. We expect that the 
loop corrections to the individual diagrams in Fig. 2 will be smaller in magnitude. However, 
two-pion exchange diagrams of chiral order u = 3 might have magnitudes comparable to our 
z/ = 2 terms since the diagrams in Fig. 2 are only part of the effective a- exchange. 

To simulate more realistic a-exchange it may also be necessary to explicitly include 
intermediate A-particles in Figs. 2 and 3, but that would require a thorough recalculation 



of many previous results. In a forthcoming paper |T^ we hope to present a detailed discussion 
of DW calculations which are required to obtain realistic cross sections for pp — ^ ppir^, as well 
as the details of a renormalization procedure relevant to an NLLO calculation. Due to large 
energy-momentum transfers involved in the NN -^ NNtt reaction, a full DW calculation 
can modify significantly the numerical results reported in this Letter. Our finding that the 
"recoil" corrections 0{m]^^) to T^^'^^ are large points to the necessity of examining the use 
of the static heavy-baryon "propagator", l/{v-p). As a first pragmatic step, one can think of 
replacing l/{v-p) with 1/ [v-p — p"^ / {2mN)) , which implies, however, an extension of HBxPT 
adopted in this work. 

We thank Shung-ichi Ando for pointing out an error in our calculation of the recoil 
correction to the pion rescattering diagram. This work is supported in part by the Na- 
tional Science Foundation, Grant No. PHYS-9602000 and by the Grant-in-Aid of Scientific 
Research, the Ministry of Education, Science and Culture, Japan, Contract No. 07640405. 
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TABLES 
TABLE I. The three sets of ci 2 3 parameters, in units of GeV~^, used in the text. 





Cl 


C2 


C3 


A 


-0.87 


3.34 


-5.25 


B 


-0.87 


4.5 


-5.25 


C 


-0.98 


3.34 


-5.25 



TABLE IL Sizes of the K type of diagrams, shown in Figs. 2-5, relative to the u = 1 pion 
rescattering diagram, Fig. 1(b). The ratio Rk defined in the text is given for three sets (A, B and 
C) of parameters (01^2,3), see the text and Table 1. 



pA pB pC 

^K ^K ^K 



K = I -0.70 -0.38 -0.53 

K = 11 6.7 3.6 5.1 

K = III -6.7 -3.6 -5.1 

K = l\ 9.5 5.1 7.2 

K = \ 0.18 0.10 0.14 

K = \l 0.14 0.08 0.11 

K = VII 2.6 1.9 2.0 
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FIGURES 



q ^ 



q 



(a) 



(b) 



FIG. 1. Diagrams contributing to pp -^ ppir^: Born term (a), the pion-rescattering term (b). 
The soUd hne denotes a nucleon (proton), the dashed hne a pion and the open elhpses denote 
initial-/final-state interactions. 
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Type I 



Type II 



Type III 



(C) 



FIG. 2. Two-pion exchange graphs of type I, II and III -the "cross-box" graphs. 
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Type IV 



FIG. 3. Pion-pion s-wave rescattering graph of type IV. 
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FIG. 4. Pion rescattering vertex corrections to Fig. 1(b) of type V and VI. 
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Type VII 



FIG. 5. The pion-rescattering graph with v =2 at rescattering vertex. 




FIG. 6. The five-point contact-interaction counter-term graph. 
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